A useful guide for gravitational wave observers to test modified gravity models 
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We present an extension of a previously suggested test of all modified theories of gravity that 
would reproduce MOND at low accelerations. In a class of models, called "dark matter emulators," 
gravitational waves and other particles couple to different metrics. This leads to a detectable time 
lag between their detection at Earth from the same source. We calculate this time lag numerically 
for any event that occurs in our galaxy up to 400 kpc, and present a graph of this possible time lag. 
This suggests that, gravitational wave observers might have to consider the possibiUty of extending 
their analysis to non-coincident gravitational and electromagnetic signals, and the graph that we 
present might be a useful guideline for this effort. 

PACS numbers: 98.80.Cq, 98.80.Hw, 04.62. +v 



I. INTRODUCTION 

The last forty years have been a very active and an 
exciting period for cosmologists. The breakthrough was 
the discovery of cosmic microwave background radiation 
(CMBR), an open window to the very early stages of 
the universe [1]. We realized that the universe is homo- 
geneous and isotropic on large scales 0, And very 
recently we learned that the universe is accelerating con- 
trary to the common behef [1, Q . 

Another very exciting observational discovery was the 
peculiar shape of the so called "galaxy rotation curve" . 
Instead of the velocity of stars orbiting spiral galaxies be- 
ing inversely proportional to the distance away from the 
center of the galaxy, it was observed to stay almost con- 
stant after certain point. Then even more suprising 
observational fact was the excessive amount of distortion 
in the background galaxies near a foreground mass 
Assuming that there is only baryonic matter in the uni- 
verse, it is not possible to explain these two observations 
with Einstein's equations, 
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Two of the following can be done to avoid this discrep- 
ancy. We can modify the right hand side by assuming 
additional dark matter, that would interact only gravi- 
tationally. Or we can modify the left hand side of the 
equation, namely the Einstein's tensor. The latter ap- 
proach suggests to modify general relativity on cosmo- 
logical scales. And it was partially achieved in 1983 by 
Milgrom 15] with the so called "MOND" for Modified 
Newtonian Dynamics. Partially, because it is not a mod- 
ification of general relativity but Newtonian dynamics on 
cosmological scales. The modification of GR would then 
be the underlying fundamental theory that would give 
MOND in the non-relativistic limit. 
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In order to understand Eq.([T]), one needs to be able 
to handle both sides of the equation. That means iden- 
tifying dark matter directly, by a non-gravitational ex- 
periment. Therefore direct detection is crucial for dark 
matter. For more than 30 years, experimentalists have 
been trying to detect dark matter without any success. 
Very recently some excitement arose about the detection, 
but turned out to be almost no signal with excluding 
new parameter space [l6|. On the other hand, the early 
relativistic [l^, and the Scalar- Tensor [l^ formulations 
of MOND led to, so called TeVeS for "Tensor- Vector- 
Scalar" , which was recently proposed by Bekenstein [l^ . 
It reproduces the MOND in the non-relativistic limit and 
does a good job on cosmological evolution, strong grav- 
itational lensing as well as post-Newtonian parameters 
(20I f2l| . But weak lensing is problematic in the case of 
colliding clusters, such as the Bullet Cluster [2^j25|. The 
observed third peak of new 5-year WMAP data[26] is also 
problematic in the context of MOND-like models. But 
still one can not claim that these models are ruled out 
(27I [28| . Thefore it is very important to find a way of 
testing these two explanations by means of observation. 



The possibility for comparing dark matter with 
MOND-like models by means of a model-independent 
test was suggested before [30l - l32| . It was proposed that 
one would see a time lag between arrival times of gravita- 
tional waves and other massless particles. And the time 
lag was calculated for three different sources; a particu- 
lar gamma-ray burst(GRB) GRB070201, the low mass x- 
ray binary (LMXB) SCo-Xl, and the supernova SN1987a. 
Other proposed tests of TeVes with gravity waves are dis- 
cussed in [29i] and references therein. In Sec. 2 we sum- 
marize the formulation of this suggested test. In Sec. 
3 we extend this calculation numerically to any source 
in Milky Way galaxy. In Sec. 4 we discuss the conse- 
quences of this effect to the externally triggered searches 
for GWs. In Sec. 5 we investigate the ambiguities in this 
calculation, and see how sensitive the time lag is to the 
choice of using different data sets. Our conclusions are 
summarized in Sec. 6. 
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II. FORMULATION 

TeVeS is an impressive achievement but a natural ques- 
tion arises when one thinks about multiple metric for- 
malisms. Why do we need to have multiple metrics? It 
was realized thatfs^ if one wants to generalize MOND, 
to get TuUy-Fisher relation[3J] and sufhcient lensing, to a 
stable, relativistic, covariant pure metric theory without 
dark matter, it would not be possible. This no-go theo- 
rem restricts the modified gravity models, that are con- 
sistent with Tully-Fisher relationship and the observed 
amount of lensing. One can violate one of the assump- 
tions of the no-go theorem by introducing another metric 
which would carry the MOND force. These class of mod- 
ified gravity models with multiple metrics, which would 
give MOND in the non-relativistic limit, were called as 
"dark matter emulators". TeVeS and the scalar- vector- 
tensor gravity theory that was proposed by Mofi^at (ssl. |36| 
are examples of dark matter emulators. These models 
all have the following property: All the ordinary matter 
should move on a geometry as if there is dark matter, 
whereas the gravity waves move on a geometry without 
any. The total time needed for a gravitational wave to 
reach the Earth from a given source is the sum of dis- 
tance divided by speed of light and additional delay due 
to gravitatiional potential of intervening matter along the 
line of sight (also known as Shapiro delay) [13, [11] . This 
can be restated as follows: Gravity waves couple to the 
usual metric of GR without dark matter and the rest of 
the matter couples to a different metric that would be 
produced with GR plus dark matter. So the answer to 
the question at the beginning of this section is, we might 
have to introduce multiple metrics for a modified grav- 
ity model in order to explain Tully-Fisher relationship 
and the observed amount of lensing consistently, without 
dark matter. 

Heuristically the calculation can be understood with a 
simple picture. Let us compare the arrival times of two 
massless particles that are produced by a distant source. 
The gravity wave is going to move on a geometry where 
there is only ordinary matter, and the other ordinary 
massless particle, which can be a photon or a neutrino, 
will move in a geometry where there is dark matter as 
well as the ordinary matter. Because of this additional 
mass the ordinary massless particle would have to travel 
a geometry as if the distance is longer, which would mean 
that gravity waves arrive earlier than all the rest of the 
massless particles. In other words,, because of this addi- 
tional mass, the ordinary massless particles feel an extra 
Shapiro delay as compared to gravity waves. Therefore 
the procedure is the following: Calculate the arrival time 
of a massless particle by solving the geodesic equation 
with a mass density being just that of dark matter, which 
would be equal to the difference of arrival times of these 
two massless particles. 

We will calculate the time lag in the context of bi- 
metric theories mimicking cold dark matter which has a 
spherically symmetric distribution. One can express the 



time lag of an event in terms of its initial (xi) and final 
{X2) positions, by solving corresponding geodesic equa- 
tions perturbatively. The result [sj for the time lag is. 
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M{r) is the mass function. 
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Therefore, in order to mimick dark matter one has to 
choose a dark matter density profile p(r) and that will 
determine the amount of time lag. 



III. CALCULATION 

In this Section, we revisit a calculation that was done 
in a previous work. The time lag of photons from 
GRB070201[32], a short duration Gamma Ray Burst 
(GRB) which is believed to have happened at the An- 
dromeda Galaxy, was calculated. To evaluate the total 
time lag, one needs to include contributions from both 
the Milky Way and Andromeda galaxies. In the men- 
tioned work, the data for the density profiles of these two 
galaxies were taken from two different groups. Therefore 
we decided to do the calculation again, based on the data 
of a group where they did a numerical work on the den- 
sity profile of both of the galaxies [1^ . 

As stated above, to calculate the time lag we need to 
know the form of the dark matter density profile. The 
NFW profile can be used for this purpose, which assumes 
the following density profile of dark matter halos in nu- 
merical simulations 14011 . 
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The parameters, that we need to calculate the time 
lag, are po the characteristic density, and ro the radius of 
the halo. The summary of the Shapiro delays for GRB 
070201 is given in Table ID One can notice that the time 
lag becomes larger with increasing virial mass of the An- 
dromeda galaxy. 

The next step is considering a source which is in the 
Milky Way and estimating the time lag by solving Eq. ^ 
for At. At this point it is useful to analyze an ap- 
proximate expression for At using the previous work of 
Woodard and Soussa [i^. Forcing the dark matter em- 
ulators to mimick dark matter isothermal halo distribu- 
tion, results into the following approximate values for 
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Data Set 


p{GeV/cm^) 




At{days) 


MW - Klypin [39] 


0.185 


1.0 X 10^^ 


426 


MW - Ascasibar [41] 


0.347 


1.0 X 10^^ 


421 


M31 - Klypin [39] 


0.188 


1.6 X 10^^ 


634 


M31 - Tempel [42] 


0.661 


1.0 X 10^^ 


383 



TABLE I: Shapiro delays for GRB 070201 from the NFW 
profiles of the Milky Way and Andromeda using two different 
data sets. 



AA{r) and AB{r): 
A{r) w 1 + e + , B{r) « 1 



e + e,ln( — ). (6) 



Here e is the usual Schwarzschild term e = {2GM / c?r) 
and the e* is the extra small parameter due to this effect, 
e, = 2v1/c^ where given the asymptotic rotation speed 
for Milky Way e, is at the order of 6 x 10~^. Since 
rs ~ 8fcpc the logarithm will be of order one. Looking 
at Eq.(l2|) and observing the direct proportionality of Ai 
with AA(r) and AB{r), one can see that the time delay 
would approximately be equal to 10^^ times the total 
time of flight in flat space. 

But our aim is to calculate the time lag numerically, 
due to the difficulty of solving the Eq. ([2]) analytically. Us- 
ing the relevant parameters for the NFW density profile 
of the Milky Way ^] , the numerical solution is depicted 
in Figure [TJ 

TimeLag 




distance(kpc) 



FIG. 1: Shapiro delays for sources located in Milky Way. 

If we look at Figure [1] we can see that the time lag for 
an event that is located at 100 kpc is approximately 100 
days, which is 10~^ times 3 x 10^ light years (the total 
time of flight of a photon in flat space) . 

The time lag would surely depend on the angular po- 
sition of the source, since we are not the center of the 
Milky Way. Therefore one would like to see the effect of 
the direction of the position of a source in the whole sky. 
This was done for a source that is located up to 400 kpc 



away from the Earth, and it is demonstrated in Figure 
[21 It can be seen that the angular position of the source 
might have a maximum of 5% effect on the Shapiro delay 
of the considered source. 



TimeLag 




FIG. 2: The angular dependence of Shapiro delays for sources 
located in Milky Way. The units of R A (Right ascension) and 
Dec(Declination) is coverted to radians. 



IV. EXTERNALLY TRIGGERED SEARCHES 
FOR GWS 

Externally triggered search of gravitational waves is a 
strategy to look for possible signals around sources which 
might produce both GWs and other particles (external 
triggers). The information provided by an external trig- 
ger is extremely useful to impose extra requirements for a 
possible GW signals as well as to increase the confidence 
of the detection. Knowing the source direction would al- 
low us to look for only relevant portion of the sky. One 
can even infer the frequency of the possible GW signal 
based on the information provided by an external trigger 
and make a data analysis at a specific band [44] . 

Gamma ray bursts (GRBs) are one of the possible 
external triggers. A neutron star(NS) black hole or a 
NS-NS merger might be a good source for a GRB. A 
highly magnetized NS might also produce short and in- 
tense gamma rays, which are called as soft gamma ray 
Repeaters (SGRs). One sixth of short GRBs are thought 
to be because of SGRs, where during S5 of LIGO, electro- 
magnetically several hundred of them are were observed 
[45| . Disruption of neutron star's crust would produce an 
increase in the rotation energy of a neutron star which 
results in increase of freqency, so called a "pulsar glitch" . 
This might also lead to emission of GWs. Last but not 
least, supernova explosions would surely produce both 
GWs and neutrinos. This triggered detection on the 
other hand, would provide us more information about 
neutrino mass as well as the mechanism of supernova ex- 
plosions. 

But still, there has not been a single direct detection 
of GWs. The search is based on the assumption that, 
the GW signal and the external trigger are coincident in 
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M31 


Mtot(lOi2j^0) 


Milky Way 


Mtot(lOi'M0) 


Corteau [46] 


1 oq+0.18 

^■•J'J-o.is 


Xue [51] 




Evans [47] 




Smith [52] 


-'-•^^-0.54 


Fardal [48] 


0.74«:1^ 


Wilkinson [53] 


1.91?;? 


Seigar [49] 


r, 70+0.02 
' 'J-0.02 


Sakamoto [54] 


J--^-0.7 


Ibata [50] 


U-'O-0.13 


Battaglia [55] 


1 tr+0.2 



TABLE II: DifTerent mass estimates of the Milky Way and 
the Andromeda with the corresponding error bars. 



time within a small time window. Therefore it is possible 
that, if "dark matter emulators" describe gravity then we 
might be looking at the wrong place. This possibility was 
suggested before and the possible time lag was calculated 
for certain events [s^j • But it is crucial for GW observers 
(especially the people who are involved with externally 
triggered search) to know how big this effect is for a par- 
ticulat object and the amount of error. In this work we 
provide an estimate of this possible time lag for any ob- 
ject located at Milky Way and current uncertainties are 
discussed in the next section. 



V. UNCERTAINTIES IN THE CALCULATION 

There are three kinds of uncertainties in the time lag 
calculation. First one is the position of the assumed 
source. The uncertainity in the angular position of the 
GRB's and other sources are determined by a relatively 
good accuracy. In a previous work For GRB-070201 the 
error was shown to be of the order of 2% [32,]. The second 
kind of uncertainty is the choice of the dark matter den- 
sity profile and this choice had an effect which wass less 
than 1% j33]- The third kind of uncertainty is the data 
input (parameters) coming from different groups who are 
working on numerical simulations of dark matter density 
profiles. And this choice has the biggest effect on the 
calculation. 

In order to calculate the time lag in Eq.® we need 
to use these parameters coming from numerical simula- 
tions. They enter into the calculation through the mass 
function. In this work, to get the mass function, we used 
NFW dark matter density profile for modified gravity 
models to mimic dark matter. One can see from Table 
m that using a data set of a different group changed the 
time lag from 383 days to 634 days for the Andromeda 
galaxy. And this difference is expected; since the virial 
mass estimation of M31 is almost twice of one another. 

Therefore the biggest uncertainty comes from the de- 
termination of the mass of the galaxies. Looking at Table 
m the results for our galaxy are not very different, unlike 
M31. But still one naturally ould start wondering about 
the case of mass estimate for our galaxy as well. Some 



of the recent mass estimates for both M31 and Milky- 
Way galaxy are shown in Table |IT1 The results of virial 
mass estimates are quite different from one another. This 
might be misleading since the virial radius which is con- 
sidered in these works are sometimes different from one 
another. But even taking that into consideration, the er- 
ror bars of these estimates are too big, varying from 13% 
to 140% ,that we can not give a precise estimate of the 
time lag for objects that are far from us. 

To summarize, our prediction of this time lag has cer- 
tain uncertainties, where the biggest one is the input 
coming from the numerical simulations of dark matter 
density profiles of the galaxies. For an event which will 
occur at a very distant location, more than lOOkpc, we 
can only conclude that the time lag will be big. Therefore 
the predictivity of this test is much better for an event 
which is closer to us lOkpc). The errors of the most 
recent mass estimates are at the order of 20% and nu- 
merical simulation community is hopeful that these mass 
estimates will get much better in recent future for our 
galaxy [sel]. 



VI. CONCLUSIONS 

A peculiar property of dark matter emulators is that 
gravitational waves couple to a different metric than the 
ordinary matter. This results in a time lag between the 
arrival times of these particles coming from the same 
source, which would be a powerful check of these models 
with a single detection of a gravitational wave. A si- 
multaneous detection of gravitational wave and a photon 
would rule out entire class of modified gravity models. 
But if a time lag is observed which is at the order magni- 
tude of the one mentioned in this work, that would mean 
general relativity is wrong and there is no dark matter. 

By a numerical analysis we made a rough estimation 
of the time lag for any source in our galaxy as far as 400 
kpc, and demonstrated it in Figure [TJ The time lag de- 
pends almost linearly on the distance of the source from 
us, and is approximately 10~^ times the total time of 
flight. Estimation of the time lag depends on the data 
input (mass function of Milky- Way) which is provided by 
numerical simulations. Therefore the accuracy of this, 
highly depends on the error bars coming from these sim- 
ulations. The current state of affairs of the numerical 
simulations of Milky- Way galaxy puts some constraints 
on our prediction. Therefore, at this stage we can con- 
clude that the best objects, for this proposed test, are 
the ones that are close to us(~ 10 kpc). 

We also have to point that Bekenstein takes the ap- 
proach as velocity of GWs is smaller than photons 19] . 
But this would lead to an immediate invalidation of 
TeVeS from the boundary of Moore and Nelsonjs^] via 
gravitational cherenkov-radiation. Therefore we do not 
agree with Bekenstein's point of view and force the "dark 
matter emulators" to agree with the observed rotation 
curve and weak lensing results. 
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Besides its own significance, gravitational wave obser- 
vation is very important since it will open new windows 
in our understanding of the universe. One can test whole 
class of alternate theories of gravity which would results 
into peculiar imprints such as additional polarization 
states. Higher order theories of gravity, Brans-Dicke the- 
ory, massive graviton theories and Chern-Simons modi- 
fied gravity are just some of the alternate theories of grav- 
ity that we can test (see [H, [11] and references therein) . 

This work should serve as a guide for gravita- 
tional wave observers who would like to search for im- 
prints of modified gravity models without dark matter. 
LIGO/VIRGO are currently doing joint science opera- 
tions at design sensitivities and searches for gravitational 
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